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Abstract 

Recent work in the literature had evaluated the energy-momentum tensor of a Casimir apparatus 
in a weak gravitational field, for an electromagnetic field subject to perfect conductor boundary 
conditions on parallel plates. The Casimir apparatus was then predicted to experience a tiny 
push in the upwards direction, and the regularized energy-momentum tensor was found to have 
a trace anomaly. The latter, unexpected property made it compelling to assess what happens 
in a simpler case. For this purpose, the present paper studies a free, real massless scalar field 
subject to homogeneous Dirichlet conditions on the parallel plates. Working to first order in the 
constant gravity acceleration, the resulting regularized and renormalized energy-momentum tensor 
is found to be covariantly conserved, while the trace anomaly vanishes if the massless scalar field is 
conformally coupled to gravity. Conformal coupling also ensures a finite Casimir energy and finite 
values of the pressure upon parallel plates. 
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I. INTRODUCTION 



Ever since Casimir discovered that suitable differences of zero-point energies of the quan- 

n 

tized electromagnetic field can be made finite and produce measurable effects (ll, several 



efforts have been produced to understand the physical implications and applications of this 
property [2j, 3l, N, \& In particular, we are here concerned with the recent theoretical dis- 

nnnn n 

covery that Casimir energy gravitates [3, |8|, 12, [10J . In Ref. |9j, some of us proved this as 
part of an investigation which led to the evaluation of the energy-momentum tensor of a 
Casimir apparatus in a weak gravitational field. In that investigation, the functional integral 
quantization of Maxwell theory was applied, with perfect conductor boundary conditions on 
parallel plates at distance a from each other. On using Fermi-Walker coordinates, where 
the (21,22) coordinates span the plates, while the z = x 3 axis coincides with the vertical 
upwards direction (so that the plates have equations z = and z = a, respectively), and 
working to first order in the constant gravity acceleration g, the spacetime metric reads as 



ds 2 = -c 2 + e^j dt 2 + dxj + dx\ + dz 2 + 0(|x| 2 ), (1.1) 

where e = ^r. The resulting regularized (with point-split method) and renormalized energy- 
momentum tensor (T^) was found to be covariantly conserved, with photon and ghost Green 
functions obeying the mixed boundary conditions of the problem and satisfying the Ward 
identities [9] . All of this would have been completely reassuring, had it not been for the fact 
that such a (T^) was found to have a trace anomaly 

W> = a" 3 / (-) , (1-2) 

where 



'G 



z\ 71 hg / z 15 cos(7rz/a) 



'1.3) 

~a) 360 c \a 2tt sm 3 (7rz/a) J v 
As is by now well known, even though the classical action is invariant under conformal 
rescalings of the metric, a trace anomaly may arise because some counterterms may occur 
which fail to possess the same invariances as the classical action (see Comments on chapter 



28 in Ref. 



111]). In four spacetime dimensions, such counterterms are quadratic in Riemann, 



Ricci and scalar curvature of the background. Moreover, if the boundary is nonempty, the 
trace anomaly receives further contributions from local invariants which are cubic in the 
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extrinsic-curvature tensor Kij of the boundary, i.e. 

K> Kj Kl K*K lm K lm , (K?f. 

All of this is made precise by the proportionality between trace anomaly and the local heat- 

leat equation for an operator of Laplace 
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kernel coefficient 02 [12( obtained by studying the 
type on a Riemannian manifold with boundary [lj 

Our analysis of the energy-momentum tensor, however, was Lorentzian rather than Eu- 
clidean. We obtained the Hadamard function as twice the imaginary part of the Feynman 
Green function, and then used the point-split method to obtain (T^). It became therefore 
important to check the consistency of the above procedure by studying a simpler problem, 
and for this purpose we have here focused on a scalar problem. Our "Casimir" appara- 
tus involves a massless scalar field in curved background, subject to homogeneous Dirichlet 
conditions on parallel plates with mutual separation a. 

Section II is a summary of basic properties of a free massless scalar field in curved space- 
time: action functional, classical energy-momentum tensor and its regularized expression. 
Section III evaluates the Feynman Green function up to first order in e. Section IV ob- 
tains the regularized and renormalized energy-momentum tensor of the quantum theory, 
while Sec. V evaluates the Casimir energy. Discussion of the results and open problems are 
presented in Sec. VI. 



II. FREE MASSLESS SCALAR FIELD IN CURVED SPACETIME 

The action functional for a free massles scalar field <fi coupled to gravity reads as 

s = -\j (<M : " + W 2 ) V^gfa, (2.1) 

where g is the determinant of the metric tensor g^, R is the scalar curvature and £ is a 
real parameter taking the value ~ in the case of conformal coupling and in the case of 
minimal coupling. On requiring stationarity of the action functional under variations of 0, 
one obtains the field equation 

(□-£R)0 = O, (2.2) 
where Q = g^V^Vj, is the wave operator. 
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The classical energy-momentum tensor is obtained from functional differentiation of the 
classical action with respect to the metric, i.e. 

2 5S 



rpl_W _ 

If use is made of the standard variational identities 



(2.3) 



5R = -R^Sg^ + g^g ua {5g pu ,, u - 8g va , pv ), 
and by introducing the anticommutator 



(x), <f){y)}+ = (j>{x)(j>{y) + 4>(y)(j>(x) 



(2.4) 



one finds eventually 



+ ^ 



,<7 ) 



+ | f rt - -g^R 



(2.5) 



At this stage, the point-split method comes into play, and we re-express every anti- 
commutator in (2.5) according to 



</>]+ = Jim i { [<^ V , 0] + + [0^, 0'] + } , 



lim [ 



(2.6) 

(2.7) 
(2.8) 



On introducing the Hadamard two-point function (hereafter, the brackets ( ) denote the 
(vacuum) expectation value) 



H(x,J) = ([<f>(x),<f>(x')] + ) ) 



(2.9) 



we get (cf. [16]) 



/rpflVK 



lim 

x'—*x 



(1-20 



W iU + H 



i - - ) <r#/ 



i<T (#./ + if./) + | (V - ~<ri?) # 



(2.10) 
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Note that the use of general equations has payed off: if we now focus on Ricci-flat space- 
times, the effect of £ remains, which would not be obvious if Ricci-flatness were imposed in 
(2.1). Hereafter, on taking the coincidence limit, we need of course the geodesic parallel dis- 
placement bivector P^, which performs parallel displacement of vectors along the geodesic 
from x' to x. In general, it is defined by the differential equations [9j 

<r»P» = ^P* = 0, (2.11) 



a(x, x') being the Ruse-Synge world function [171 ]. equal to half the geodesic distance between 
x and x', jointly with the coincidence limit 



lim Pt = 



K- (2-12) 



Equation (2.11) means that the covariant derivatives of P^ ; vanish in the directions tangent 
to the geodesic joining x and x'. Thus, the bivector P^,, when acting on a vector B v at x', 
gives the vector P M , which is obtained by parallel transport of B v> to x along the geodesic 
connecting x and x', i.e. 

B 1 * = P», B v ' . (2.13) 
In our problem, we have from Eq. (1.1) the metric tensor 

<^ = diag(-l-e-,l,l,l), (2.14) 

V (X * 

with contravariant form 

gT ~ diag (-1 + s-, 1, 1, l) + 0(e 2 ). (2.15) 

The orthonormal tetrad e a ji relating the spacetime metric g^ v to the Minkowski metric r] a b 
according to 

9r = e a / uVab (2.16) 

is thus given by 

e\ = e\ = 8 1 , V* = 1, 2, 3. (2.17) 

The bivector P^, hence reads as [l^l 

K = 9" P W b p e a u , ~ diag (l + |-(*' - z), 1, 1, l) + 0(e 2 ), (2.18) 

while 

P/ = 9^P% ~ diag (l + ±-{z - z'), 1, 1, l) + 0(e 2 ). (2.19) 



The corresponding geodesic is described by the equations 

t = constant, x\, x^ = constant, 
z^ = if)(z' - z) + z, if) e [0, 1]. 



(2.20) 
(2.21) 



III. FEYNMAN GREEN FUNCTION TO ZEROTH AND FIRST ORDER 

The equations II were rather general, whereas now we consider, following the introduction 
and our previous work [h]] , Fermi- Walker coordinates for a system of two parallel plates 
in a weak gravitational field. To first order in the e parameter of Sees. I and II, the only 
nonvanishing Christoffel symbols associated with this metric are therefore 

+ (3., 

We are now in a position to evaluate the scalar counterpart of the analysis in Ref. [9J, i.e. 
we compute the wave operator the Feynman Green function of the hyperbolic operator 
( | | — £R), and eventually the Hadamard function and the regularized energy-momentum 
tensor. 

Indeed, a Green function of the operator ruling the field equation (2.2) obeys the differ- 
ential equation 

{U-iR)G(x,x') = - 5 -^£X. (3.2) 

The Feynman Green function Gp is the unique symmetric complex- valued Green function 
which obeys the relation [3] 

5G = G5F G, 

where F is the invertible operator obtained from variation of the action functional with 
respect to the field. This definition is well suited for the purpose of defining the Feynman 



Green function even when asymptotic flatness does not necessarily hold [191 ] . 

In our first-order expansion in the e parameter, the scalar curvature gives vanishing 
contribution to Eq. (3.4), which therefore takes the form (hereafter Q° = rj^d^d^) 

□° + + ^rfir) *') = (3.3) 



{a + ez)d£ w (a + ez)dzJ v ' ' y/=g 
We now follow our work in Ref. [9j and assume that the Feynman Green function admits 
the asymptotic expansion 

G F (x, x') ~ G (0) (x, x') + eG {1 \x, x 1 ) + 0(e 2 ). (3.4) 



It is a nontrivial property that this asymptotics should hold at small e for all x,x', no 
matter how close or distant from each other are the two points. Its existence is proved by 
the calculations described hereafter. Indeed, by insertion of (3.4) into (3.3) we therefore 
obtain, picking out terms of zeroth and first order in e, the pair of differential equations 

U°G (p) (x,a?) = J ( ®(x,x'), (3.5) 

\J°G^(x,x') = J^(x,x'), (3.6) 

having set 

J^(x,x') = -5(x,x'), (3.7) 

jmM ^± JM .^ + L^ GmM . (3. 8) 

Our boundary conditions are Dirichlet in the spatial variable z. Since the full Feynman 
function Gf{x,x') is required to vanish at z = 0,a, this implies the following homogeneous 
Dirichlet conditions on the zeroth and first-order terms: 



G (0) (i,s') 



= 0, (3.9) 

2=0, a 



G^(x,x')\ =0. (3.10) 

1 2=0, a 

To solve Eqs. (3.5) and (3.6), we perform a Fourier analysis of and G^\ which remains 
meaningful in a weak gravitational field by virtue of translation invariance. In such an 
analysis we separate the z variable, i.e. we write (cf. j^]) 

G^{x,x') = ! ^!^ 7 (o) ( ^^) e ^-(^-^)-^°(-o-^) ; (3 n) 

J (27r) d 

and similarly for G^(x,x'), with a "reduced Green function" j^-\z,z') in the integrand 
as a counterpart of the zeroth-order Green function j^(z,z') in (3.11). Equations (3.5) 
and (3.6) lead therefore to the following equations for reduced Green functions (hereafter 



A V^fci): 

J^ + A 2 ) 7 ( °W) = -<W), (3-12) 

J + *) = + g« - ~) 7<W)- (3.13) 

By virtue of the Dirichlet conditions (3.11), reads as 

(0 )( z>) = _sin(A* < )sin(A(* > -a)) 
Asin(Aa) 



where z < = mxn(z,z'), z > = max(z,z'). The evaluation of the reduced Green function 
is slightly more involved. For this purpose, we distinguish the cases z < z' and z > z', and 
find the two equations 



+ A a ) 7 £W) = j£ 1 W), 



where 



dz 1 



^ = _lAcos(A,)-2^sin(A,) ^ _ . f ^ < ^ 



2a 



A sin(Aa) 



(i) 1 A cos(A(2 — a)) — 2zk% sin(A(z — a)) . . 

3+ — o , " sm(Az ) if z > z . 



(3.15) 



(3.16) 



(3.17) 



2a Asin(Aa) 
We have therefore two different solutions in the intervals z < z' and z > z' . In this case the 
differential equation (3.13) is solved by imposing the matching condition 



t! 1 W) = 7?V, 



jointly with the jump condition 



dz 



7+ 



d (i) 
7- 



<9z 



2' 

2a' 



(3.18) 



(3.19) 



Equation (3.18) is just the continuity requirement of the reduced Green function j^'(z, z') 
at z — z', while Eq. (3.19) can be obtained by integrating Eq. (3.13) in a neighborhood of 



z , since 



lim -^-7 ( - 1 ^ 
*o dz 



z'+e 



lim 



y+e z z' 
5(z, z')dz 



(3.20) 



Bearing in mind Eq. (3.14) we can therefore write, for all z, z' 

1 



7 W 



4aA 2 



(kt-\>)(z + z')-kZ(z^ + z'^ 



9_ , .,2 9 

dz 



7 ( °W) 



, 9 9 sin(Az) sin(Az') 



sin (Aa) 



(3.21) 



IV. REGULARIZED AND RENORMALIZED ENERGY-MOMENTUM TENSOR 



In the previous section we have focused on the Feynman Green function Gp because it is 
then possible to develop a recursive scheme for the evaluation of its asymptotic expansion 
at small e. However, as is clear from Eq. (2.10), we eventually need the Hadamard function 
H(x,x'), which is obtained as 9J 



H(x,x') = 2lmG F (x,x') ~ 2Im(G (0) (x, x') + eG (1 \x, x')) + 0{e 2 ). 



(4.1) 



The coincidence limits in (2.10), with the help of (2.13) and (2.18), make it necessary to 
perform the replacements 



H^ v + H,, vl -> P/H. yu + P/H.^, H/ -> g^P/H, ap/ , H,^ -> P/p/h.,^, (4.2) 



Hence we get the asymptotic expansion at small e of the regularized energy-momentum 
tensor according to (hereafter we evaluate its covariant, rather than contravariant, form) 



(^)~(TW)+e(T«) + 0(e 2 ), 
where, on defining s = nz/a, s' = irz'/a, we find 



and 



(T<$) = 



71 



lim 



7T 



1440a 4 s'— *s 



3-2 sin 2 s 
sin 4 s 







A 








o\ 









-1 








1 










-1 









\o 










/l 








o\ 







-1 


















-1 





1 


\o 













oo / 



7T 



1440a 4 sin 4 s 



311 637 1, 

— 7r — s H (43tt - 81s) cos 2s 

40 40 10 v ; 



s — 3tt 

+ — — — cos 4s + 5 sin 2s + 2{ix — s)s(sin 2s — 6 cot s) 



40 



+ u- 



7T 



[2(tt + s)(2 



cos 2s) H — sin 2s 
; 2 



6 J 48a 4 sin 4 s 

]7TS 
— lim — — — , 
s'^s 2a 4 (s — s ) 4 



(4.3) 



(4.4) 



(4.5) 



<2ff> = 



7T 



7T - 2S + 



sin 2 s 



(^2(tt - 2s) 



7200a 4 

+ cots (5 + 2(tt - s)s - 6(tt - s )— 1-)) 
V sin s// 



2 + 



sin 2 s 



+ U-i 



7T 



sin 

11 (7T — s)s — 1) COS S 



6 J 96a 4 sin 5 s 
+ ((n — s)s + 1) cos3s — 2(7r — 2s)(3sins + sin3s)j , 



(4.6) 
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(Iff) 



7T 



7TS 



7T COS S 



(4.7) 
(4.8) 



1440a 4 ' 720a 4 ' V" 5 16a 4 sin 3 s' 

The next step of our analysis is the renormalization of the regularized energy-momentum 

tensor. For this purpose, following our work in Ref . , we subtract the energy-momentum 

tensor evaluated in the absence of bounding plates, i.e. 

/l o\ 



lim 



7T 



2a 4 (s 



-A4 



and 



lim 



TVS 



2a 4 (s - s') 4 



0-100 
0-10 
\0 3/ 

A o o o\ 



\0 0/ 



(4.9) 



(4.10) 



To test consistency of our results we should now check whether our regularized and renormal- 
ized energy-momentum tensor is covariantly conserved, since otherwise we would be outside 
the realm of quantum field theory in curved spacetime, which would be unacceptable. In- 
deed, the condition 

W M (T I »,) = Q (4.11) 
yields, working up to first order in e, the pair of equations 



^(T 3 ( 3 0) ) = 0, (e°tenn) 



which are found to hold identically for all values of £ in our problem. 

The trace of (T M1/ ) is obtained as (g^ u being the contravariant metric in (2.15)) 



(4.12) 
(4.13) 



from which we find a ^-dependent part 

3tt 2 (2 + cos2s) 



CI 



0(e< 



(4.14) 



7T 



la 4 sin 4 s 



32a 4 sin 5 s 

— (1 + (-7T — s)s) cos 3s + 2(7r — 2s) (3 sins + sin 3s 



1 — 11(71" — S)s) COS S 



(4.15) 
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Interestingly, the value £ = | which yields conformal invariance of the classical action (2.1) 
with wave equation (2.2) is the same as the value of £ yielding no trace anomaly, unlike what 
happens for Maxwell theory in Ref. 9J], where the conformally invariant action is found to 
lead to a trace anomaly in quantum theory with mixed boundary conditions. 



V. CASIMIR ENERGY AND PRESSURE 



In order to evaluate the energy density p of our "scalar" Casimir apparatus, we project 
the regularized and renormalized energy-momentum tensor along a unit timelike vector 



0,0,0 . This yields 



V~9oo ' 

p = (T^u 
-2 + - 



7T 



+ 



7T 



1440a 4 7200a 4 



— 37T + 6s + 



10 



sin 2 s 



2(tt - 2s) 



x 



sin 2 s 



+ cots (5 + 2(vr - s)s + 6 



S(-7T + S, 

sin 2 s 



1\ f vr 2 (2 + cos2s) 
' " ~ 6/i 8a 4 sin 4 s 



7T 



192a 4 sin 5 s 



— 5 + 22(7r — s)s ) coss 



+ ^5 + 2(7r — s)sj cos 3s — 4(7r — 2s) (3 sins + sin 3s) 
The energy .E stored within our Casimir cavity is given by 



E 



d 3 £, 



-gp, 



(5.1) 



(5.2) 



where d 3 S is the volume element of an observer with four-velocity u^, and V c is the volume 
of the cavity. The integration in (5.2) requires the use of approximating domains, i.e. the 
^-integration is performed in the interval ((,a — (), corresponding to -(£, a — Q in the s 
variable, taking eventually the ( — > limit. We thus obtain 



E f 



7l 2 A 



ir Ae 



1440a 3 5760a 3 



■kA / e\ cosC 
^~6j4a?\ + V C™sin 3 C' 



(5.3) 



where A is the area of parallel plates. Note that the conformal coupling value £ 



is 



picked out as the only value of £ for which the Casimir energy remains finite. In this case, 
reintroducing the constants h, c and writing explicitly e, we find 



Tick 2 A ( 1 ga 
1440 a 3 V + 2^ 



In the same way, the pressure on the parallel plates is found to be 9j 

^ 2 n 2 e ( 1 



= o) 



7T 



+ 



480a 4 1440a 4 



7re cos s 
lim 



6 7 16 a 4 "-^6 sin 3 s ' 



(5.4) 



(5.5) 
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r, / \ vr 2 7r 2 e ( . l\ TIE ,. coss „. 

P £ ( z = a) = + + U - - j lim — 5-. (5.6) 

CV ; 480a 4 1440a 4 V 6/ 16a 4 sin 3 s V ; 

Once again, one can get rid of divergent terms by setting £ = i which leads to 



TV 2 he f 2 qa \ 



n 2 he ( ^ 2 ga \ . . , 

p ^ = °) = -«o7 I 1 -3?)- (5 ' 8) 



VI. CONCLUDING REMARKS 



In our paper we have obtained, for the first time in the literature, a detailed evaluation of 
the regularized and renormalized energy-momentum tensor for a scalar Casimir apparatus in 
a weak gravitational field, the plates being parallel plates upon which a real massless scalar 
field is required to obey homogeneous Dirichlet conditions. Moreover, the trace anomaly 
is found to vanish provided the free scalar field is conformally coupled to gravity. The 
conformal coupling also ensures a finite Casimir energy and finite pressure on the plates. 

The electromagnetic case is definitely more involved: the boundary conditions are a mix- 
ture of Dirichlet and Robin conditions, and if the Casimir apparatus is set in a weak gravi- 
tational field the trace anomaly is not found to vanish [9j , despite that the classical Maxwell 
action is conformally invariant in four spacetime dimensions. At least three investigations 
are now in order: 

(i) To repeat the scalar analysis with Robin boundary conditions, to understand whether 
the latter are responsible for the nonvanishing trace anomaly found in Ref. for the 
electromagnetic case. 

(ii) To work out the relation (if any) between our small-e asymptotics of the Feynman Green 
function GV, with the following limit and the usual approach where one first 
considers the Schwinger-DeWitt asymptotics of Gf at small values of the world function 
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20 



2l|. The two approaches are not obviously equivalent nor easily comparable, but 
the novel features found in Ref. |&] make it compelling to produce further efforts along both 
lines. 



(iii) To relate our results to the enlightening energy-momentum analysis in Ref. 22] . 
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